A fast - Monte Carlo toolkit on GPU for treatment plan dose recalculation in proton therapy

In the context of the particle therapy a crucial role is played by Treatment Planning Systems (TPSs), tools aimed to compute and optimize the tratment plan. Nowadays one of the major issues related to the TPS in particle therapy is the large CPU time needed. We developed a software toolkit (FRED) for reducing dose recalculation time by exploiting Graphics Processing Units (GPU) hardware. Thanks to their high parallelization capability, GPUs significantly reduce the computation time, up to factor 100 respect to a standard CPU running software. The transport of proton beams in the patient is accurately described through Monte Carlo methods. Physical processes reproduced are: Multiple Coulomb Scattering, energy straggling and nuclear interactions of protons with the main nuclei composing the biological tissues. FRED toolkit does not rely on the water equivalent translation of tissues, but exploits the Computed Tomography anatomical information by reconstructing and simulating the atomic composition of each crossed tissue. FRED can be used as an efficient tool for dose recalculation, on the day of the treatment. In fact it can provide in about one minute on standard hardware the dose map obtained combining the treatment plan, earlier computed by the TPS, and the current patient anatomic arrangement

Novel imaging and quality assurance techniques for ion beam therapy: a Monte Carlo study

Ion beams exhibit a finite and well defined range in matter together with an “inverted” depth-dose profile, the so-called Bragg peak. These favourable physical properties may enable superior tumour-dose conformality for high precision radiation therapy. On the other hand, they introduce the issue of sensitivity to range uncertainties in ion beam therapy. Although these uncertainties are typically taken into account when planning the treatment, correct delivery of the intended ion beam range has to be assured to prevent undesired underdosage of the tumour or overdosage of critical structures outside the target volume. Therefore, it is necessary to define dedicated Quality Assurance procedures to enable in-vivo range verification before or during therapeutic irradiation. For these purposes, Monte Carlo transport codes are very useful tools to support the development of novel imaging modalities for ion beam therapy. In the present work, we present calculations performed with the FLUKA Monte Carlo code and preliminary experimental studies

The FLUKA code for application of Monte Carlo methods to promote high precision ion beam therapy

Monte Carlo (MC) methods are increasingly being utilized to support several aspects of commissioning and clinical operation of ion beam therapy facilities. In this contribution two emerging areas of MC applications are outlined. The value of MC modeling to promote accurate treatment planning is addressed via examples of application of the FLUKA code to proton and carbon ion therapy at the Heidelberg Ion Beam Therapy Center in Heidelberg, Germany, and at the Proton Therapy Center of Massachusetts General Hospital (MGH) Boston, USA. These include generation of basic data for input into the treatment planning system (TPS) and validation of the TPS analytical pencil-beam dose computations. Moreover, we review the implementation of PET/CT (Positron-Emission-Tomography / Computed- Tomography) imaging for in-vivo verification of proton therapy at MGH. Here, MC is used to calculate irradiation-induced positron-emitter production in tissue for comparison with the +-activity measurement in order to infer indirect information on the actual dose delivery

Diamond detectors for dose and instantaneous dose-rate measurements for ultra-high dose-rate scanned helium ion beams

Background: The possible emergence of the FLASH effect—the sparing of normal tissue while maintaining tumor control—after irradiations at dose-rates exceeding several tens of Gy per second, has recently spurred a surge of studies attempting to characterize and rationalize the phenomenon. Investigating and reporting the dose and instantaneous dose-rate of ultra-high dose-rate (UHDR) particle radiotherapy beams is crucial for understanding and assessing the FLASH effect, towards pre-clinical application and quality assurance programs. Purpose: The purpose of the present work is to investigate a novel diamond-based detector system for dose and instantaneous dose-rate measurements in UHDR particle beams. Methods: Two types of diamond detectors, a microDiamond (PTW 60019) and a diamond detector prototype specifically designed for operation in UHDR beams (flashDiamond), and two different readout electronic chains, were investigated for absorbed dose and instantaneous dose-rate measurements. The detectors were irradiated with a helium beam of 145.7 MeV/u under conventional and UHDR delivery. Dose-rate delivery records by the monitoring ionization chamber and diamond detectors were studied for single spot irradiations. Dose linearity at 5 cm depth and in-depth dose response from 2 to 16 cm were investigated for both measurement chains and both detectors in a water tank. Measurements with cylindrical and plane-parallel ionization chambers as well as Monte-Carlo simulations were performed for comparisons. Results: Diamond detectors allowed for recording the temporal structure of the beam, in good agreement with the one obtained by the monitoring ionization chamber. A better time resolution of the order of few μs was observed as compared to the approximately 50 μs of the monitoring ionization chamber. Both diamonds detectors show an excellent linearity response in both delivery modalities. Dose values derived by integrating the measured instantaneous dose-rates are in very good agreement with the ones obtained by the standard electrometer readings. Bragg peak curves confirmed the consistency of the charge measurements by the two systems. Conclusions: The proposed novel dosimetric system allows for a detailed investigation of the temporal evolution of UHDR beams. As a result, reliable and accurate determinations of dose and instantaneous dose-rate are possible, both required for a comprehensive characterization of UHDR beams and relevant for FLASH effect assessment in clinical treatments

A fast - Monte Carlo toolkit on GPU for treatment plan dose recalculation in proton therapy

In the context of the particle therapy a crucial role is played by Treatment Planning Systems (TPSs), tools aimed to compute and optimize the tratment plan. Nowadays one of the major issues related to the TPS in particle therapy is the large CPU time needed. We developed a software toolkit (FRED) for reducing dose recalculation time by exploiting Graphics Processing Units (GPU) hardware. Thanks to their high parallelization capability, GPUs significantly reduce the computation time, up to factor 100 respect to a standard CPU running software. The transport of proton beams in the patient is accurately described through Monte Carlo methods. Physical processes reproduced are: Multiple Coulomb Scattering, energy straggling and nuclear interactions of protons with the main nuclei composing the biological tissues. FRED toolkit does not rely on the water equivalent translation of tissues, but exploits the Computed Tomography anatomical information by reconstructing and simulating the atomic composition of each crossed tissue. FRED can be used as an efficient tool for dose recalculation, on the day of the treatment. In fact it can provide in about one minute on standard hardware the dose map obtained combining the treatment plan, earlier computed by the TPS, and the current patient anatomic arrangement

The FLUKA Code: Developments and Challenges for High Energy and Medical Applications

n/

The role of nuclear reactions in Monte Carlo calculations of absorbed and biological effective dose distributions in hadron therapy

Monte Carlo codes are rapidly spreading among hadron therapy community due to their sophisticated nuclear/electromagnetic models which allow an improved description of the complex mixed radiation field produced by nuclear reactions in therapeutic irradiation. In this contribution results obtained with the Monte Carlo code FLUKA are presented focusing on the production of secondary fragments in carbon ion interaction with water and on CT-based calculations of absorbed and biological effective dose for typical clinical situations. The results of the simulations are compared with the available experimental data and with the predictions of the GSI analytical treatment planning code TRiP

Development and benchmarking of a dose rate engine for raster-scanned FLASH helium ions

Background: Radiotherapy with charged particles at high dose and ultra-high dose rate (uHDR) is a promising technique to further increase the therapeutic index of patient treatments. Dose rate is a key quantity to predict the so-called FLASH effect at uHDR settings. However, recent works introduced varying calculation models to report dose rate, which is susceptible to the delivery method, scanning path (in active beam delivery) and beam intensity. Purpose: This work introduces an analytical dose rate calculation engine for raster scanned charged particle beams that is able to predict dose rate from the irradiation plan and recorded beam intensity. The importance of standardized dose rate calculation methods is explored here. Methods: Dose is obtained with an analytical pencil beam algorithm, using pre-calculated databases for integrated depth dose distributions and lateral penumbra. Dose rate is then calculated by combining dose information with the respective particle fluence (i.e., time information) using three dose-rate-calculation models (mean, instantaneous, and threshold-based). Dose rate predictions for all three models are compared to uHDR helium ion beam (145.7 MeV/u, range in water of approximatively 14.6 cm) measurements performed at the Heidelberg Ion Beam Therapy Center (HIT) with a diamond-detector prototype. Three scanning patterns (scanned or snake-like) and four field sizes are used to investigate the dose rate differences. Results: Dose rate measurements were in good agreement with in-silico generated distributions using the here introduced engine. Relative differences in dose rate were below 10% for varying depths in water, from 2.3 to 14.8 cm, as well as laterally in a near Bragg peak area. In the entrance channel of the helium ion beam, dose rates were predicted within 7% on average for varying irradiated field sizes and scanning patterns. Large differences in absolute dose rate values were observed for varying calculation methods. For raster-scanned irradiations, the deviation between mean and threshold-based dose rate at the investigated point was found to increase with the field size up to 63% for a 10 mm × 10 mm field, while no significant differences were observed for snake-like scanning paths. Conclusions: This work introduces the first dose rate calculation engine benchmarked to instantaneous dose rate, enabling dose rate predictions for physical and biophysical experiments. Dose rate is greatly affected by varying particle fluence, scanning path, and calculation method, highlighting the need for a consensus among the FLASH community on how to calculate and report dose rate in the future. The here introduced engine could help provide the necessary details for the analysis of the sparing effect and uHDR conditions

On the parametrization of lateral dose profiles in proton radiation therapy.

Abstract Purpose The accurate evaluation of the lateral dose profile is an important issue in the field of proton radiation therapy. The beam spread, due to Multiple Coulomb Scattering (MCS), is described by the Moliere's theory. To take into account also the contribution of nuclear interactions, modern Treatment Planning Systems (TPSs) generally approximate the dose profiles by a sum of Gaussian functions. In this paper we have compared different parametrizations for the lateral dose profile of protons in water for therapeutical energies: the goal is to improve the performances of the actual treatment planning. Methods We have simulated typical dose profiles at the CNAO (Centro Nazionale di Adroterapia Oncologica) beamline with the FLUKA code and validated them with data taken at CNAO considering different energies and depths. We then performed best fits of the lateral dose profiles for different functions using ROOT and MINUIT. Results The accuracy of the best fits was analyzed by evaluating the reduced χ2, the number of free parameters of the functions and the calculation time. The best results were obtained with the triple Gaussian and double Gaussian Lorentz–Cauchy functions which have 6 parameters, but good results were also obtained with the so called Gauss–Rutherford function which has only 4 parameters. Conclusions The comparison of the studied functions with accurate and validated Monte Carlo calculations and with experimental data from CNAO lead us to propose an original parametrization, the Gauss–Rutherford function, to describe the lateral dose profiles of proton beams

A new PET prototype for proton therapy: comparison of data and Monte Carlo simulations

Ion beam therapy is a valuable method for the treatment of deep-seated and radio-resistant tumors thanks to the favorable depth-dose distribution characterized by the Bragg peak. Hadrontherapy facilities take advantage of the specific ion range, resulting in a highly conformal dose in the target volume, while the dose in critical organs is reduced as compared to photon therapy. The necessity to monitor the delivery precision, i.e. the ion range, is unquestionable, thus different approaches have been investigated, such as the detection of prompt photons or annihilation photons of positron emitter nuclei created during the therapeutic treatment. Based on the measurement of the induced β+ activity, our group has developed various in-beam PET prototypes: the one under test is composed by two planar detector heads, each one consisting of four modules with a total active area of 10 × 10 cm2. A single detector module is made of a LYSO crystal matrix coupled to a position sensitive photomultiplier and is read-out by dedicated frontend electronics. A preliminary data taking was performed at the Italian National Centre for Oncological Hadron Therapy (CNAO, Pavia), using proton beams in the energy range of 93–112 MeV impinging on a plastic phantom. The measured activity profiles are presented and compared with the simulated ones based on the Monte Carlo FLUKA package